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We develop a temperature dependent empirical pseudopotential theory to study the electronic 
and optical properties of self-assembled quantum dots (QDs) at finite temperature. The theory 
takes the effects of both lattice expansion and lattice vibration into account. We apply the theory 
to the InAs/GaAs QDs. For the unstrained InAs/GaAs heterostructure, the conduction band offset 
increases whereas the valence band offset decreases with increasing of the temperature, and there 
is a type-1 to type-11 transition at approximately 135 K. Yet, for InAs/GaAs QDs, the holes are 
still localized in the QDs even at room temperature, because the large lattice mismatch between 
InAs and GaAs greatly enhances the valence band offset. The single particle energy levels in the 
QDs show strong temperature dependence due to the change of confinement potentials. Because of 
the changes of the band offsets, the electron wave functions confined in QDs increase by about 1 - 
5%, whereas the hole wave functions decrease by about 30 - 40% when the temperature increases 
from to 300 K. The calculated recombination energies of exciton, biexciton and charged excitons 
show red shifts with increasing of the temperature, which are in excellent agreement with available 
experimental data. 

PACS numbers: 68.65.Hb, 73.22.-f, 78.67.Hc 
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I. INTRODUCTION 

During the past two decades, enormous progress has 
been achieved in understanding the electronic and opti- 
cal properties of self-assembled quantum dots (QDs) both 
through theory and experiments, stimulated by their po- 
tential applications in QDs laser at room temperature^i"— 
and as qubits and quantum photon emitters at low 
temperature 4^— For the former applications, it gener- 
ally requires high density and highly uniform QDs. The 
QD laser has been demonstrated with much lower thresh- 
old current Jc and much higher material and differential 
gains as compared to the semiconductor quantum well 
lasersi^ For the latter applications, the preparation of 
single QD is crucial. A number of methods have demon- 
strated the feasibility to isolate single QD from QDs 
ensemble fiSrJ^ Rabi oscillation of exciton and charged 
excito n^^i^^ in single QD have been observed experimen- 
tally, showing that the charge and spin quantum states 
in single QD can be coherently controlled via optical 
method. The single and entangled photon emission from 
single QD have also been demonstrated experimentally^^ 
that are much brighter— than the traditional parameter- 
down entangled photon source.— i^^ These experimental 
achievements pave the way for future application of QDs 
in quantum computation. 

On the other hand, the development of the atom- 
istic theories, including the empirical pseudopotential 
method^S"— and the tight-binding models^^"— provide 
deep insight to the electronic and optical properties of 
self-assembled QDs. The atomistic theories of QDs not 
only give results that agree well with experiments^^"— 
but also greatly improve our understanding of the prop- 



erties of QDs. The atomistic models capture the correct 
point group symmetry of the QDs, which is missing in 
the continuum model. Therefore they can give correct 
interpretation of some subtle properties of the QDs, e.g., 
the fine structure splitting (FSS)i^ and light polariza- 
tion of excition^i Unfortunately, so far all the theories 
of QDs have been restricted to zero temperature. 

Temperature is a very important degree of freedom 
in experiments to tune the electronic and optical prop- 
erties of QDs. For example, in QDs laser, the tem- 
perature has be used to tune the laser wavelengthii"— 
In the QD-cavity system, the temperature is gener- 
ally used to tune the resonance between the QDs and 
the cavity in order to achieve strong coupling between 
the two quantum systems i^"— The temperature depen- 
dent optical spectra of (single and ensemble) QDs have 
been investigated intensively in experiments in the past 
decadeSfi2ri^iiiii^i^i2£"iS New physics, for example, the 
formation of excitonic polaron22i^d£ which is due to 
strong coupling between exciton and optical phonons, 
may be found in QDs at high temperature. However, 
a theoretical understanding of the temperature effects in 
QDs is still missing. Therefore, to facilitate the future 
device applications of QDs, the development of a temper- 
ature dependent theory is not only of theoretical interest, 
but also of practical importance. 

In this work, we develop such a temperature depen- 
dent atomistic pseudopotential theory to study the elec- 
tronic and optical properties of QDs at finite tempera- 
ture. We take the effects of both lattice expansion and 
lattice vibration into account. The latter is done by 
introducing a temperature dependent dynamical Debye- 
Waller factor to the structure factor. We first examine 



the temperature dependent electronic structures of bulk 
In As and GaAs, and then apply the theory to investigate 
the electronic and optical properties of self-assembled 
InAs/GaAs QDs. The calculated temperature dependent 
photoluminescence (PL) spectra of QDs are in excellent 
agreement with available experimental data. 

The rest of the paper is organized as follows. In Sec. [ll] 
we introduce the temperature dependent empirical pscu- 
dopotential method (TDEPM). In Sec. |llll we study the 
electronic structures of bulk InAs, GaAs using TDEPM, 
including the energy band gaps and band offsets, etc. 
We present the temperature dependent band offsets for 
InAs/GaAs QDs in Sec. IIV Al and the single particle 
energy levels and wave functions of InAs/GaAs QDs in 
Sec. IIVBI We discuss the temperature dependent optical 
spectra of InAs/GaAs QDs in Sec. llVCi and summarize 
in Sec. M 
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FIG. 1: (color online) Experimental temperature dependent 
lattice constants for InAs and GaAs. The data are taken from 
Ref. H 



II. METHODOLOGY 

We consider InGaAs QDs embedded in the center of 
a 60 X 60 X 60 GaAs 8-atom unit cell. Periodic bound- 
ary condition is used to obtain the single particle energy 
levels. To study the electronic and optical properties of 
the QDs at finite temperature, we introduce temperature 
dependent pseudopotcntials in the single-particle Hamil- 
tonian, 
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Va{v -R„a,e,T), 
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where Va{i^, e, T) is the strain dependent screened empir- 
ical pseudopotential for atom of type a and atom index 
n at temperature T. R„q is the optimized atom position 
from valance force field (VFF) method.— i^ In a strained 
lattice, the atomic potential is assume to have the form 
of, 



w„(r,e,T) = w„(r,T)[l + 7„Tr(e(r))] 



(2) 



where Tr(e(r)) is the local hydrostatic strain at r. 7^ is 
fitted to the deformation potentials of the bulk materials. 
The atomistic theory naturally captures the correct point 
group symmetry of the QDs even at high temperature. 

The effects of lattice vibration can be taken into ac- 
count using Brooks- Yu theory^^ which has been applied 
to study the energy gap of bulk materials in the context 
of pseudopotentials4^"— The total atomic potential at 
temperature T, which is the sum of all atomic potentials 
Va{Y,e,T), can be written as, 

V{v,T) =(^^«„(q)5(q)e"i-)T 

q a 



q a 



where q is the reciprocal lattice vector. Va(c\} is the 
Fourier transform of the screened atomistic potential at 



zero temperature, which takes the form of;^ 

ao{q^ - ai) 



Va{q) 



Q,2e"39=' - 1 



(4) 



where q ~ |q|, and the parameters ao, ai, q;2, as are 
fitted to the bulk properties of InAs and GaAs, includ- 
ing the band gaps, band offsets, effective masses, etc. 
{Sa{<\))T is the averaged structure factor over all the 
phonon configurations at temperature T , 



{SM))T = {Y.e-'^^-^--)T. 



(5) 



Assuming R„q = R^j„ + u^ , where u^ is the amplitude 
of the phonon mode. For any q, we have 



/e-iq-"-)j,^e-3((q-"-)'>T. 



(6) 



using Wick's theorem.— Therefore, the temperature ef- 
fects to the atomic potentials are equivalent to consider 
a temperature dependent structure factor. 
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in Eq. [31 where Wa (q, T) is the dynamical Debye- Waller 
factor for the a-th clement. 



2I^„(q,r) = ((q.u„)2)j.. 



(8) 



For simplicity, we assume the system to be isotropic, then 
we have 



Ty„(q,r) = i|qp(^,2) 



(9) 



where (w^) is the total mean-square displacement for 
atom of type a at temperature T, including the contri- 
bution from acoustic (A) and optical (O) phonons. 



("a) = ("a)A+ (Wa)o- 
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Wc use Dcbyc model for acoustic phonons 
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where 04 = 9h/Muj]j, and S/j = /icj^i is the Debye 
temperature. The second term in the bracket is from 
zero point quantum fluctuation, which can be absorbed 
to the zero-temperature pseudopotentials. At high tem- 
perature, the acoustic phonon displacement {v?)a oc T, 
whereas in the low-temperature limit, {v?)a oc T^. For 
the optical phonons, the average displacement can be 
written as 



{u\ 



Oih 



^ru^o/i<B 



1 



(12) 



where wq is the frequency of the optical phonon. We ne- 
glect the dispersion of optical phonons here. Generally, 
the contribution of the optical phonon could be signif- 
icant only at high temperature.— Although 04 and as 
can be calculated directly using the parameters of bulk 
materials, the approximations made during the deriva- 
tions of Eq. pT|) and ([T2|) can introduce some errors 
to the energy gap of semiconductors at finite tempera- 
ture. To overcome this problem, we treat Q!4 and a^ as 
fitting parameters, which are fitted to the temperature 
dependent band structures in combination with the zero- 
temperature empirical pseudopotentials. Therefore, the 
dynamical Debye- Waller factors might be different from 
the real physical Debye- Waller factors of the system. 



TABLE I: Band parameters obtained from the empirical 
pseudopotential calculations compared with the experimental 
values.— 1^ Energies are in unit of eV. m%, mj!^^[100], m^,j[lll] 
and m*i^ [100] are effective masses of electron and holes at F 
point, flgap and apis^ denote the deformation potential of the 
band gap and the Fi5„ point, h is the valance band biaxial de- 
formation potential. Ao and Ai are the spin-orbit splittings 
at the Fi5„ and I/i„ points respectivity. 





GaAs 


InAs 


Property 


EPM 


Expt. 


EPM 


Expt. 


-C/gap 


1.528 


1.52 


0.423 


0.42 


Ex,^ 


-2.763 


-2.96 


-2.365 


-2.40 


Ex,^ 


1.937 


1.98 


2.069 


2.34 


Ex,^ 


2.232 


2.50 


2.514 


2.54 


Els^ 


-1.041 


-1.30 


-0.872 


-1.26 


El,^ 


2.232 


1.81 


1.568 


1.71 


ml 


0.067 


0.067 


0.023 


0.023 


mlulWO] 


0.316 


0.40 


0.371 


0.35 


mL[lll] 


0.825 


0.57 


0.986 


0.85 


m,\[100] 


0.092 


0.082 


0.029 


0.026 


dgap 


-7.879 


-8.33 


-6.804 


-5.7 


"ri5„ 


-1.110 


-1.0 


-0.829 


-1.0 


b 


-1.567 


-1.7 


-1.631 


-1.7 


Ao 


0.362 


0.34 


0.384 


0.39 


Ai 


0.201 


0.22 


0.286 


0.27 



TABLE IL Fitted pseudopotential parameters for InAs/GaAs 
in Eq.(121) and Eq.Q. A plane- wave cutoff of 5 Ryd is used. 



parameters 


Ga 


As (GaAs) 


In 


As (InAs) 


ao 


476845.70 


11.9753 


771.3695 


26.8882 


ai 


1.9102 


3.0181 


1.6443 


2.9716 


02 


22909.50 


1.1098 


18.1342 


1.2437 


as 


0.1900 


0.2453 


0.3940 


0.4276 


Jc 


2.5215 


0.0 


2.1531 


0.0 


Ctso 


0.1035 


0.0976 


0.5973 


0.0976 



TABLE III: Debye temperatures,-^ the optical phonon ener- 
gies and the fitted 04, as parameters for GaAs and InAs. 





eD(K) 


cation 


anion 


Bulk 


aJo(meV) 04 as 


ct4 as 


GaAs 
InAs 


344 
247 


35.36 0.3024 0.0786 
29.6 0.1014 0.0984 


0.1530 0.0024 
0.0828 0.0084 



To determine the temperature dependent pseudopo- 
tential, we first determine the pseodupotential parame- 
ters ao - C(3 and 7 at zero temperature by fitting them 
to the electronic structures of bulk materials, including 
the effective mass, and energies of the high symmetry F, 
X and L points, etc. The target values and fitted values 
for GaAs and InAs are compared in Table HI which are 
in good agreement. The parameters ap - a^ and 7 are 
presented in TablelTTl With these parameters at hand, we 
then determine the values of a4 and a^ by fitting them 
to the temperature dependent energy gap of bulk materi- 
als, which can be well described by the empirical Varshni 
formulai^2i^ 



AEg{T) = - 



T + C2' 



(13) 



where ci and C2 are the Varshni parameters. For GaAs, 
ci = 0.5405 meV/K and C2= 204 K and for InAs, ci = 
0.276 meV/K and C2 = 93 KM The fitted parameters for 
a4 and as are summarized in Table. IIIIl 

The single-particle Hamiltonian ([Ij can be solved by 
expanding the wave functions into a linear combination 
of Bloch bands (LCBB)^ 



^^=E 



'n,k,AV'n,e,k,A,T, 



(14) 



.,k,A 



where ipn.eM,\.T is the bulk Bloch bands with orbital n 
and wave vector k close to F point at finite temperature 
T, and A = (InAs, GaAs). The experimental lattice con- 
stants for InAs and GaAs at given temperature, as shown 
in Fig. [U are given as input to construct the Bloch ba- 
sis. At each temperature, we relax the dot-|-matrix struc- 
ture using VFF method to get the the atomic position 

Due to the spatial confinement, the carries in the QDs 
have strong Coulomb interactions. The many-particle 



Hamiltonian read as, 



H 



Yl e^^lV"* 



ijkl 



(15) 



where tpi = CiV'i(r) is the field operator with correspond- 
ing single particle energy e^. F^,- are the Coulomb inte- 
grals, 



ij " 



drdr' 



e(r — r')|r — r' 



(16) 



Here, e(r — r') is the screened dielectric functioni^ The 
many-particle Hamiltonian is solved using a configura- 
tion interaction method)^ where the many-particle wave 
functions are expanded on the Slater determinants con- 
structed from the confined electron and hole levels. This 
method has been successfully applied to studying the 
electronic and optical properties of InAs/GaAs QDs and 
the obtained results are in very good agreement with the 
experimental observations j^^j^Sii^i^ 
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FIG. 2: (color online) The band structures of GaAs at T = 
(a) and T = 300 K (b) calculated by the TDEPM. 



III. TDEPM FOR BULK MATERIALS 

We first test our method for bulk materials. Figure [21 
depicts typical band structures of GaAs at T = and T 
= 300 K. The overall band structures are quite similar 
to those at zero temperature even at rather high temper- 
ature (300 K). However, the energies of high symmetry 
fc-points F, A, and L have different response to the tem- 
perature. In principle, all the energies of these fe-points 
should be taken as the target values to determine the val- 
ues of a4 and a^ . Unfortunately, the experimental data 
of the energies of these fc-points at finite temperature 
are not available, therefore, tentatively, we fit the poten- 
tials only to the temperature dependent energy gaps at F 
point. The potentials can be improved by fitting to the 
energies of more fc-points in the future. 

We present in Fig. [3] the change of energy gap AEg 
as a function of temperature for GaAs. The filled circles 
are the results taken from experiments;^ whereas the 
solid line is calculated from TDEPM. The fitting error is 
less than 0.2 meV in the whole temperature range. The 
dashed line is the change of energy gap of GaAs taken 
account only lattice vibration, whereas the dotted line 
is the result with only lattice expansion. As we see, to 
accurately describe the red shift of the band gap with 
respect to temperature, both the lattice expansion and 
lattice vibration have to be taken into account in the the- 
ory. These results agree with the results obtained from 
temperature dependent tight-binding method by Pour et 
aliS Similar features are also found for InAs using the 
parameters given in Table IIIIl 

The conduction band offset {AEc) and valance band 
offset (AEy) between InAs and GaAs, defined as 
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FIG. 3: (color online) The red shift of band gap as a func- 
tion of temperature for GaAs. The filled circles represent the 
experimental data. The dashed line is the result calculated 
from TDEPM with only lattice vibration whereas the dotted 
line is the result with only lattice expansion. The solid line 
is from TDEPM taking account of both lattice vibration and 
lattice expansion. 



are very important to the electronic structures of the 
InAs/GaAs hetcrostructures, because they are important 
for the confinement of electron and hole in InAs/GaAs 
QDs. We show the change of CBM and VBM of InAs 
and GaAs with temperature in Fig. |4] (a), and the tem- 
perature dependent band offsets between InAs and GaAs 
in Fig. 11(b). We find: 

(i) The CBM and VBM of GaAs change much larger 
than their counterparts in InAs. 

(ii) The CBM generally changes much larger than 
VBM. At high temperature (T > 100 K), CBM and VBM 
change approximately linearly with respect to tempera- 
ture. The changes of the band offsets with respect to the 
temperature are presented in Fig. |4](b) for the VBM and 
in the inset for CBM. For electron, the temperature effect 



generally enhances the band offset, whereas for hole, the 
band offset decreases with the increasing of the tempera- 
ture. Interestingly, there is a type-I to type-II transition 
at T = 135 K for unstrained InAs/GaAs. However, be- 
cauase of the strain effects, the holes are still localized in 
InAs/GaAs QDs (see Sec. lIVp . In what follows we will 
show that the change of band offsets at finite temperature 
will greatly modify the electronic and optical properties 
of QDs. 
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FIG. 4: (color online) (a) Changes of CBM and VBM as 
functions of temperature for InAs and GaAs bulk materials, 
(b) The valence band offset and the conduction band offset 
(inset) between InAs and GaAs as functions of temperature. 



IV. TDEPM FOR INAS/GAAS QDS 

In Sec. mil we study the temperature dependent elec- 
tronic structures of bulk InAs and GaAs. The band off- 
sets between InAs and GaAs are greatly modified due to 
the temperature effects and there is a type-I to type-II 
transition in the unstrained InAs/GaAs heterostructure. 
The change of band offsets will significantly change the 
corresponding electronic and optical properties in QDs. 
However, the temperature dependent properties are more 
complicated in QDs because of the strain effects. In this 
section, we investigate the temperature dependent elec- 
tronic and optical properties of InAs/GaAs QDS using 
TDEPM. We study QDs with different sizes and alloy 



TABLE IV: The alloy composition and size of the lens shaped 
lui^Gai-ajAs/GaAs QDs used in this work. 



# 



Base (nm) 



Height (nm) 



QD-A 


1.0 


25 


3.5 


QD-B 


0.6 


25 


3.5 


QD-C 


0.7 


25 


3.5 


QD-D 


0.8 


25 


3.5 


QD-E 


1.0 


25 


5.5 


QD-F 


0.6 


25 


5.5 


QD-G 


0.7 


25 


5.5 


QD-H 


0.8 


25 


5.5 



compositions. The alloy compositions of selected lens- 
shaped QDs are presented in Table IIVI In most of the 
cases, we use QD-A to illustrate the main physics. The 
results of other QDs will also be presented for compari- 
son. 



A. Temperature dependent band offset 
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FIG. 5: (color online) (a) Typical band profiles in strained 
QDs. HH, LH and SO correspond to heavy hole, light hole 
and spin-orbit bands respectively, (b) The band offsets at 
the center of InAs/GaAs QDs as functions of temperature 
calculated from TDEPM. 

We first investigate the strained band offsets in 
InAs/GaAs QDs, which is crucial for the electronic and 
optical properties of QDs. In previous works,— i^ the 
Bir-Pikus model is used to obtain the strain modified 
band profiles. However, the temperature dependent pa- 
rameters for the Bir-Pikus model is generally unavailable. 
Therefore, we calculate the band profiles directly using 
the TDEPM. After the lattice relaxations for the dot sys- 



tern, we construct the 8-atom unit cell according to the 
local strain, and then calculate the band structures us- 
ing the TDEPM. Typical band offsets in InAs/GaAs dots 
along the [100] direction at T = 150 K is compared to 
those of zero temperature in Fig. [5] (a). We see that the 
overall profiles of the heave hole (HH) and light hole (LH) 
and spin-orbit (SO) bands arc still quite similar at the 
two temperatures. The SO band is lower than the HH 
and LH by about 400 meV in the matrix and is greatly 
enhanced in the dot materials. The degeneracy of the HH 
and LH bands is broken because of the biaxial straiui ^^i^" 
The strained band offsets of CBM and VBM in the cen- 
ter of InAs/GaAs QDs are presented in Fig. [5] (b) in 
the temperature range of - 300 K. For electrons, the 
confinement is enhanced with increasing of the tempera- 
ture, whereas for holes, the confinement decreases from 
320 to 180 mcV. However, unlike the bulk materials, even 
at high temperature, AE^, is always positive, indicating 
that the hole is always localized in the QDs. The change 
of the band offsets greatly modifies the electronic and 
optical properties of QDs, as shall be discussed below. 



Temperature dependent single particle levels 
and wave functions 
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FIG. 6: (color online) Single-particle energy levels of the con- 
fined states as functions of temperature in QD-A. The shaded 
regime at the top (bottom) panel corresponding to the un- 
confined states with energies higher (lower) than the CBM 
(VBM) of GaAs. 

The evolution of single particle energy levels as func- 



tions of temperature is illustrated in Fig. |6l We show the 
results for QD-A here. Similar features are also found for 
all other QDs. We show all confined electron states, and 
the highest 20 hole states. Because of the enhancement 
of confinement potential for electrons, more states are 
confined in the dots as the increasing of the tempera- 
ture. The labels of e^ (hi) represent the energy levels 
of electron (hole) in ascending (descending) order. One 
can also use angular momentum s, p, d, etc. to label the 
wave functions. For instance, the states cq and ho can be 
labeled as s and the ei_2 (/ii,2) states are usually labeled 
by pi,2, etc. 
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FIG. 7: (color online) The single-particle s-p level energy 
spacing as functions of temperature for electrons and holes 
in the InAs/GaAs QDs. 

The s-p energy level spacing is shown in Fig. [71 which 
is defined as. 



SCsp = 



•^pi 



(18) 



for both electrons and holes. Because of the increase 
of confinement potential, the electron s-p level spacing 
increases with the increasing of the temperature. In con- 
trast, for holes, the energy level spacings decreases with 
the increasing of the temperature due to the decrease of 
the hole confinement potential. For QD-A, we find that 
the level spacing of electron increases from 64 to 70 meV 
when the temperature increases from K to 300 K. For 
hole, the energy difference decreases from 11 to about 1 
meV. The change of level spacings in QDs may be mea- 
sured from PL emission spectra. 

The p-orbit splitting is another important quantity for 
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FIG. 8: (color online) The single-particle p-orbit energy split- 
tings as a functions of temperature for electrons and holes in 
InAs/GaAs QDs. 



the single particle levels, which is defined as, 



5epp = |e 



pi 



^V2\ 



If the QDs have C4^ or T^j symmetry, the p-orbit splitting 
is exactly zero. However, for real QDs, the highest sym- 
metry is C2„, and in alloyed QDs, the symmetry is Ci, 
the splitting is nonzero. The results of the p-orbit split- 
ting are presented in Fig. [5] for different types of QDs. 
For QD-A and QD-E, the electron p-orbit splittings are 
almost independent of temperature, whereas for other 
two types of QDs the p-orbit splittings slightly increase 
with the increasing of temperature. The results for holes 
are very different from those for electrons, as shown in 
Fig. [8](b). Because the hole level spacings are very small 
and the anti-crossing between the hole levels may occur 
when increasing the temperature, the p-level splittings 
are not monotonic functions of the temperature. At high 
temperature, the p-level spacing may even exceed the 
s-p level spacing, which will never happen for electrons. 
The p-orbit splitting can be measured experimentally via 
pump-probe spectroscopy.— i^ 

The change of confinement potentials may also change 
the shapes of the wave functions. We present the squared 
envelope wave functions of electrons and holes at T = 0, 
100 and 300 K in Fig. [9] for QD-A. The number on the 
lower right corner of each small panel represents the per- 
centage of the density confined in the QDs. For eo,i,2, 
we find that the confined densities increase by about 1% 
- 2%, whereas for the higher states, the confined densi- 
ties may increase about 4% - 5%. Although the electron 
squared wave functions confined in QDs slightly increase. 



T = OK 


T=100K 


T = 300 K 


eO 


hO 


eO 


hO 


eO hO 


'»; 


■m 


■m 


m 


* v" 


82.3% 


88.8% 


82.7% 


86.7% 


83.8% 52.4% 


el 


hi 


el 


hi 


el hi 




o 


V 

w 


o 




78.8% 


81.6% 


79.5% 


78.5% 


81 .2% 51 .9°/. 


e2 


h2 


e2 


h2 


e2 h2 




o 


V 


O 


W 


79.0% 


77.4% 


79.7% 


71 .8% 


81 .4% 56.0°/. 


e3 


h3 


e3 


hS 


e3 h3 


- V 


o 




^ 


v^ 


75.3% 


82.2% 


76.3% 


80.9% 


78.8% 41 A% 


e4 


h4 


e4 


h4 


e4 h4 










•> c> 


75.2% 


70.0% 


76.2% 


65.8% 


78.7% 38.6°/, 


e5 


h5 


e5 


h5 


e5 h5 


- 


^ 


- 


' 




75.0% 


76.8% 


76.2% 


69.6% 


79.0% 44.7°X 



FIG. 9: (color online) The squared envelope wave functions of 
the confined electron and hole states at T = K, 100 K and 
300 K. The number on the lower right corner of each small 
panel represents the percentage of the state densities confined 
in the QDs. 



their overall shapes hardly change. In contrast, for holes, 
the confined densities reduce dramatically when the tem- 
perature increases. For instance, the confined density de- 
creases from 88.8% to 52.4% for h^ when the temperature 
increases from to 300 K. For /15, the confined density 
decreases from 76.8% to 44.7%. The shapes of the en- 
velope functions for holes also change dramatically. For 
instance, the Iiq state is Gaussian-like at low temper- 
ature, but at high temperature there is a node at the 
center of the wave function. This change of the wave 
function is due to the enhancement of the interfacial ef- 
fect in the QDs, because the confinement is very small 
for holes at high temperature. At low temperature only 
the tall QDs have such interfacial hole statesi ^^'^^ The 
change of the wave functions with respect to temperature 
can be measured experimentally using magnetotunneling 



spectroscopy] 
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C. Temperature dependent PL emission spectrum 

In this section, we investigate the temperature depen- 
dence of the optical spectra of InAs/GaAs QDs. The op- 
tical spectra of single InAs/GaAs QDiSi^i^iiiii^i^ have 
been measured at low temperature. The highest mea- 
suring temperature as far as we know is up to 100 K, 
performed by Ortner et alw^ At higher temperature, the 
signal to noise ratio may become very low and the single 
QD emission is usually hard to detect. However, for QDs 
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FIG. 11; (color online) The red shifts of emission lines of (a) 
X, (b) XX, and (c) X" and (d) X^ as functions of tempera- 
ture. The symbols are the results calculated from TDEPM, 
whereas the solid lines are the fitted results using Varshni 
formula. 



FIG. 10: Temperature dependent PL emission spectra of QD- 
A. The emission lines are broaden by 0.5 meV with a Lorentz 
function. 



ensemble, the emissions from QDs can be resolved even 
at room temperature. Therefore in this work, we study 
the optical spectra of QDs up to 300 K. 

Figure [TU] depicts the PL emission spectra of QD-A 
from K to 300 K. The energy difference between S and 
P shell emission can be approximated by the sum of s-p 
single particle energy level spacing of electron and hole, 
i.e. 



AE 



SP 



5el 






Although the change of level spacing in Fig. [7] can not 
been directly measured from PL emission spectra, the 
sum of them can be measured. The energy difference 
between S and P shells for QD-A is 80 meV at zero 
temperature and reduces to 72.5 meV at 300 K. 

Figure [TT] depicts the temperature dependent energies 
of the primary exciton, biexiciton and charged excitons 
for dots A, C, E, G. All the energies of exciton complexes 
show red shift as the temperature increases. For all dots, 
we find that the red shifts of the exciton complexes emis- 
sion lines can be fitted very well using Varshni formula, 
with the fitting errors generally less than 1 meV. The 
fitting errors are slightly larger than that for the bulk 
materials (< 0.2 mcV), but much smaller than the total 
red shift of the emission line which is about 100 meV. 
This suggests that the red shift of the emission energies 
is proportional to T^ at low temperature. The values 
of the Varshni parameters are summarized in Table |Vl 
These parameters are very different than those of bulk 
In As, GaAs, and also vary from dot to dot. 

The temperature dependent optical spectra has been 



studied experimentally by several groupsi^^^ For exam- 
ple, red shift of exciton emission line in single QDs have 
been investigated by Ortner et a^ from T = K to T= 
100 K. In Fig. [I2](a), we compare our theoretical results 
(QD-B, QD-F) with the available experimental datai2& 
for the Ino.6Gao.4As/GaAs QDs. The red shift of the 
emission line of exciton agree well with QD-B (less than 
1 meV). In the inset of Fig. [121 we present the exciton 
energies as a functions of temperature, which also show 
excellent agreement. The red shift of exciton energy in 
Ref. |3a can be well described by Varshni formula using 
ci = 0.4419 meV/K and C2 = 221.77 K, with error less 
than 1 meV. This is also in a good agreement with the 
theoretical values for QD-B given in Table [V] 

In Fig. [12] (b) , we compare the theoretical results of 
QD-A and QD-C with the experimental results for QDs 
ensemble measured by Yeo et ali^ where the peak energy 
of S'-shell is chosen as the exciton emission lines at each 
temperature. For QDs ensemble, the S'-shell can be well 
resolved even at room temperature. We see that the red 
shifts of QDs ensemble A2 and A3 agree well with the 
theoretical prediction of QD-A. Moreover, the exciton 
energy of A2 and A3 at zero temperature is around 1.02 
- 1.03 cV, also agree well with the exciton energy of QD-A 
listed in Table [V] 

Wc note that some experiments^ suggest that the red- 
shift of the exciton energies proportional to y^ ^66,67 jjj_ 
stead of T^ as suggested by TDEPM at very low temper- 
ature [T < 10 K). The discrepancy may come from two 
reasons. First, it is because we do not have high accu- 
rate temperature dependent band gaps to fit at very low 
temperature at present stage. It may also partly comes 
from the approximations wc made in deriving the tem- 
perature dependent pscudopotentials. Nevertheless, in 
this temperature range, the change of exciton energies is 



TABLE V: Fitted Varshni 
parameters, Ci is in unit of 

Hvl 



parameters for different QDs. Eo is the energy at zero temperature in unit of eV. The Varshni 
meV/K and C2 is in unit of K. The geometry and alloy composition of the QDs are listed in Table 
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FIG. 12: (color online) (a) Comparison of the red shifts of 
the exciton emission energies of QD-B and QD-F to the ex- 
perimental data of single QDs (dotl, dot2) in Ref. [3a . (b) 
Comparison of the red shifts of the exciton emission energies 
of QD-A and QD-C to the experimental data of QDs ensem- 
ble {A2, A3) in Ref l3a . In both figures, the insets depict the 
exciton energies as functions of temperature. 



very small, and the difference between the experimental 
values and the theory is very subtle. 

We also calculate the FSS and the polarization of the 
mono-exciton at finite temperatures. We find that the 
FSS and the polarization are generally insensitive to the 
temperature. For example, the change of FSS is usually 



less than 1 /ieV, and the change of polarization is less 
than 5 degree when increasing the temperature from to 
100 K. This result suggests that the FSS can not been 
tuned using temperature effect. 



V. SUMMARY AND CONCLUSIONS 

We develop a temperature dependent empirical pseu- 
dopotential theory, and apply it to study the electronic 
and optical properties of self-assembled InAs/GaAs 
quantum dots (QDs) at finite temperature. The the- 
ory takes the effects of both lattice expansion and lattice 
vibration into account. The pseudopotentials correctly 
reproduce the temperature dependent band gap of bulk 
III-V semiconductors such as InAs, and GaAs, etc. We 
find that for the unstrained InAs/GaAs heterostructure, 
the conduction band offset increases whereas the valence 
band offsets decreases with the increasing of the tempera- 
ture, and there is a type-I to type-II transition at approx- 
imately 135 K. Yet, for InAs/GaAs QDs, the holes are 
still localized in the QDs even at room temperature be- 
cause the large lattice mismatch between InAs and GaAs 
greatly enhances the valence band offset. The single par- 
ticle energy levels in the QDs show strong temperature 
dependence due to the change of confinement potentials. 
As a consequence, more electron states are confined at 
higher temperature. Because of the changes of the band 
offsets, the electron wave functions confined in QDs in- 
crease by about 1 - 5%, whereas the hole wave functions 
decrease by about 30 - 40% when the temperature in- 
crease from to 300 K. The calculated recombination 
energies of exciton, biexciton and charged excitons show 
red shift with the increasing of the temperature, which 
are in excellent agreement with available experimental 
data. We expect the theory can facilitate the future de- 
vice applications of QDs. 
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